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thickness in X^^ model 



Domain wall - type solution with oscillating thickness in a real, 
scalar field model is investigated with the help of a polynomial approx- 
imation. We propose a simple extension of the polynomial approxima- 
tion method. In this approach we calculate higher order corrections 
to the planar domain wall solution, find that the domain wall with 
oscillating thickness radiates, and compute dumping of oscillations of 
the domain wall. 
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Abstract 



1 Introduction 



Topological defects constitues an important class of solutions in field-theoretical 
models with degenerate vacua. They play very important role in several 
branches of physics. Let us mention here field-theoretical cosmology and the 
cosmic strings hypothesis (see 0, 0), dynamics of superconductors, su- 
perfiuids and liquid crystals in condensed matter physics (see [|| - 0) as well 
as a flux tube in QCD. 

This short and incomplete list shows the necessity of having effective com- 
puting methods to study the dynamics of topological defects. In spite of 
the increasing development of mathematical techniques to solve nonlinear 
equations, exact solutions seem not to be the rule and numerical methods 
have been the most common approach to study properties of topological de- 
fects. Therefore analytical and perturbative methods are of great interest 
and importance. 

In this paper we study domain wall type solution in the model. In 
paper excited domain wall of this kind was considered and the radia- 
tion emitted from the domain wall was found. In our approach we further 
develop the method of a so-called polynomial approximation (see 0, ||T0|), 
which is used to construct the domain wall type solution with time-dependent 
thickness. We show how to compute corrections to the polynomial solution. 
The corrections consist of two parts: the static one and the time-dependent 
one. The time-dependent part of the correction contains radiation emitted 
by domain wall with oscillating thickness. 

The plan of our paper is the following. In the next section we present 
the model. In section ^ we derive the time- dependent planar domain wall 
solution in the polynomial approximation. Section |^ is devoted to the de- 
tailed analysis of this solution. We present the method of finding correction 
to the polynomial solution. In sections ^ and ^ we calculate static and time- 
dependent part of this correction. In section |^ we analyze the backreaction of 
the radiation emitted on the domain wall. In section ^ we shortly summarize 
the main points of our work. 
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2 The Model 



We consider the model with single, scalar, real - valued field defined by 
the action: 

S = J d^x[-^v,ud^<l>d''<^ 

where 77^,^ = diag{—l, 1, 1, 1) and A, v are positive constants. The corre- 
sponding equation of motion for the field $ has the form: 

df.d^'^ - 2A$(<I>2 - v^) = 0. (2) 

The energy functional for the model is given by the formula: 

= d^xid^^d"^ + A($2 - v^]. (3) 

It is convenient to rescale the space - time coordinates and the scalar field as 
follows: 

V 

t = ax^, 
^ = ax^, 



where a = ykiP. The new variables are dimensionless. The vacuum values 
of in the considered model are equal to ±1. Configuration of the field 
which smoothly interpolates between these two vacua is called the domain 
wall. Our goal is to construct the domain wall configuration, localised on 
the x^ — x"^ plane, with the time-dependent width. The planar domain wall 
distinguishes the direction perpendicular to the wall plane, given in our case 
by the coordinate lines of C,- As we consider the configuration with (p inde- 
pendent of the coordinates x^ and x^, we can restrict our approach to the 
1 + 1 dimensional model. Equation (Q) then takes the form: 

Static solutions of this equation are well-known. They have the form: 

0(O = ±tanhe, (6) 
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3 Domain wall solution in the polynomial ap- 
proximation. 

In this section we construct approximate domain wall type solution of Eq. 

with time- dependent width. In order to realize this we use the method 
of a polynomial approximation, whose detailed description is given in p. 
The basic idea of this approach is to approximate the scalar field inside 
the domain wall by the polynomial in the variable ^ with time-dependent 
coefficients. Thus, inside an interval [—^1,^0] (^0 and ^1 are positive) we 
have: 

The domain wall solution is characterized by the fact, that for sufficiently 
large | C, \ the field approaches its vacuum values ±1: 

= +1 for e>eo, 

= -1 for e<-^i- (8) 

It is possible (see 0) to adopt more accurate asymptotics for large | ,^ |, with 
exponential correction exp(— 2q;i^) to the vacuum values ±1. For simplicity 
of further calculations we use expressions (H). 

We can tune accuracy of our approximation changing degree of the polyno- 
mial (^. One can easily notice, that the cubic Ansatz presented above is 
the simplest, nontrivial choice. The Ansatz (0) should be smoothly matched 
with the vacuum solutions at ^ = ^0 and ^ = —^1. The matching conditions 
follow in a standard manner from Eq. (|^). One integrates Eq. over ^ in 
arbitrary small intervals [—^1 — e, —^1 + e] ; [^0 ~ Cj^o + e] and lets e — > 0. 
This implies: 



From Eq. (H) one gets: 



0(t,eo) = +1, 
= -1. 



(9) 
(10) 



The Ansatz (^ with matching conditions (P) - ([T0|) allows us to find solu- 
tion of Eq. d^) in the proposed form. Inserting the expansion (|^) into the 
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matching conditions we get, analogously to the case of cylindrical domain 
wall discussed in 0, the following conditions for functions a(t), b(t), c(t) and 
parameters Ci- 

^0 = ^1, a{t) = ^, b{t) = 0, c(t) = -^. (11) 

Then we insert the expansion (J^) into @ and equate to zero coefficients 
in front of successive powers of ^. After some easy algebra, considering 
conditions (|lT]) we get the following equation: 

d + -a^ -2a = 0. (12) 
9 ^ ^ 

Function c(t) is related to a(t) (see formula ([TT|) ) and the last step in our 
approach is to solve Eq. (1121). We construct an approximate solution of Eq. 



(12), which is convenient for the further analysis. With variables redefined 



as follows: 

a{t) = ^Ait), T = 2t, (13) 

Eq. ( [T^ ) takes the form: 

A = A-A^, (14) 

where now dot denotes the derivative with respect to r. Numerical analysis 
of Eq. (|l^ shows that there exist periodic solutions of two kinds. Periodic 
solutions of the first kind oscillate around static solution A = 1 and take only 
positive values; periodic solutions of the second kind oscillate around static 
solution A = —1 and take only negative values. We can easily analyse the 
oscillating solutions of Eq. ( p!4D rewriting it in the form of the two equations of 
the first order by substitution A = B. Then periodic solutions are equivalent 
to closed trajectories in (A,B) configuration space. (1,0), (—1,0) are the 
central critical points. The (A,B) configuration space is presented in Fig. 1. 
In the further part of our discussion we consider only the periodic solutions 
oscillating around A = 1, because of the connection between A and the 
domain wall thickness parameter, ^o- When A{t) > for all r we have also 
.^0,^1 > 0, which is consistent with the interpretation of these parameters. 
The periodic solutions oscillating around A = —1 are symmetric (in the sense 
that when one changes and by each other it has no infiuence on the 
dynamics of our system), so we don't discuss this situation. Considering all 
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these remarks we restrict our investigation to solutions oscillating with small 
amplitude around A = 1. We can write: 



A{t) = 1 + ^{t), (15) 
where ^I/(r) is small, periodic function. Inserting (|l^) into Eq. (0) one gets: 



^ + ^ 



(16) 



Expression (|T6D has the form of oscillator equation with nonlinear terms. Our 
goal is to find its perturbative solution with initial condtions: ^1/(0) = \&0) 
\&(0) = 0, which describes initially static, sqeezed wall. One can solve the 
nonlinear oscillator equation of type ( ]T6| ) using Krylov - Bogolubov method 
(see [jl9|).The general solution has the form: 



^(r) = ucosip{T) H — ao{u!) + ^ dki^) cos kip^r) + Pk{^^) smkip{T), (17) 
2 k>i 

where functions u, ip, a, (3 can be calculated in a standard manner (see 

+ 0{5'), (18) 



1^) with initial condition ti;(0) = 5. In our case we find: 



5 cos 



52 

. + -r cos 
4 4 



IH r 

16 / 



where 5 is a small parameter. Thus we get solution a{t): 



a{t) 



1 + 5 cos 



2 1 + 



16 



52 

— + — cos 
4 4 



4 1 + 
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+0(5^). 
(19) 

Formula (|19[) agrees well with the numerical solutions of Eq. (p!2D. 
In the last step we have to insert the solution ([T9|) into the Ansatz (0) con- 
sidering conditions ([11]). The final expression for the field of the oscillating 
domain wall in the polynomial approximation has the form: 



i - + -5 cos nt (i - i^) + COS 2VLt U - 3f 



L2 



©(eo(t)- 1 ^ I) 



^m-Ut))-Q{-Ut)-0 + o{5' 



(20) 
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where 

n = 2 + ^5^ + 0(6''). (21) 

o 

The function ^o(^) can be regarded as the half- width of the domain wall in 
the ^ coordinate. It is given by the formula: 

^o(t) = + + ^6\os2nt + 0{S^). (22) 

From Eq. (|2T|) one can get period of oscillation of the domain wall T = 2tt/Q 
: T = IT — 0(5^). In the case of linear oscillations around static solution (one 
can regard this situation neglecting nonlinear terms in (0) - it is admissible 
because is small by definition) one gets Q = 2 and T = 2. The formula 
( pOD can be adopted to the special, static case of domain wall solution, taking 
6 = 0. It is equivalent to the (1,0) central, critical point of Eq. (|T^) . The 
general solution (p^^^t,^) can be then treated as a small oscillation (with 
amplitude given by the parameter 6) around the static solution 0^°''(O • 

#(O = ^(^-^e')0(l-|el). (23) 

Accordingly to Eq. (^OD initially the domain wall is a bit sqeezed. Then it 
oscillates with small amplitude around the static solution. 

4 Correction to the polynomial solution 

In this section we propose a simple extension of the pure polynomial solution 
obtained above. Formula (^) gives us the simplest, approximaty domain 
wall - type solution in our model. One can study wider range of phenomena 
(e.g. nontrivial asymptotics of the domain wall solution as well as radiation, 
which can be emitted by the oscillating domain wall) considering higher - 
order corrections to this solution. Our goal is to calculate perturbatively 
corrections to the zeroth order polynomial solution. Let us denote it by 
Thus, we can write the field satisfying Eq. (|]) as: 

cPit,0 = <P^'\t,0 + <P^'Kt,0- (24) 

In order to find the correction (j)^^\t,^) let us insert expression ( p^ into the 
field equation (^. Neglecting nonlinear terms in the field (f)^^^ we get: 

- 2[30(°)2 _ i]0(i) = -5^5^0(0) + 2[0(°)2 _ i]0(o)_ ^25) 
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We expect that global character of the domain wall solution during time 
evolution won't change because of topological stability (in the other words: 
the solution remains the planar domain wall type all the time, only small 
corrections can occur). That's why the linearization in (^) may be done. 
Eq. (|25| ) is a linear one with source term j{t,^) given by the formula: 



+ 2[0(°)2_i]0u 



(26) 



As it was done for the polynomial solution 
and time - dependent part: 



we can split j(t, i^to statical 



(27) 



From Eq. 

JsiO = 



27 

T 



one gets: 



0(l-U|)-(5'(e-l) + '5'(e + l)), (28) 



where prim denotes derivative with respect to ^, and: 



--6n\^ - cos fit - 3(5(4^ - ^^) cos fit 



81 



+ tH^"3^'' (e- a cos fit 



+ 

+ 



oim- 1 e I) 
mut)- 1 e I) - e(i- 1 e 



5fiM ^ - cos nt -36sgn^{l- 3^^) cos fit 



+ lsgnai-e) mo{t)-\^\)-S{l-\a)) 

+ ^5(e-acosfit 6'{ut)-\a) 

+5fi2 COS fit m + ut)) - - urn 

- [5'{^ - Ut)) - b\i - 1)] - \b\i + i,{t)) - b\i + 1)] . 



(29) 
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Eq. (|25| ) is complicated due to unpleasant form of the polynomial solution 
given by Eq. (^). We can simplify the operator on the left-hand side of Eq. 
( P^D noticing, that (p^^^ can be split into statical and time-dependent part, 
and predominant component tanh ^ (which is exact, statical solution of the 
initial field equation) can be extracted from the static part as follows: 

0W=tanhe + 5#(O + 50?^(i,O, (30) 

Inserting expression (|30|) into the evolution equation and rewriting some 
terms on the right-hand side we get finally: 

- 2 [3 tanh^ e - 1] = jsiO + Jd{t, 
+6f0(°)2_tanh2^) 0« 



+6 {6(pf'> + 2 + 2 tanh ^) 0^^) . (31) 

The component in front of (p^^^ in the third term on the right - hand side is 
nonzero only on the border of the domain wall and quickly tends to zero when 
I ^ |— >■ 00. Analogously, the component in front of (j)^^^ in the fourth term 
is of order of the small parameter 6. Thus, we have the following, strongly 
suggested method of solving Eq. (^): in the first step we put (p^^^ = on 
the right-hand side of Eq. ( pT]) and solve linear equation for 0^^^ with source 
term jsiC) + 0- The solution obtained in this way can be inserted to the 
right-hand side of Eq. (0), so in the second step we have to solve equation 
with the same linear operator but a new source term. We find full solution as 
a result of such as iterating procedure. Nevertheless if components in front 
of (p^^^ on the right-hand side are small (and that is our case), we can find 
predominant part of full solution in the first step. 

Due to the form of source term (which consist of static and time-dependent 
part) it seems natural to split (p^^^ into two components: the static one and 
the time-dependent one. They are denoted respectively by 0^^^ and ^^^^ 



From Eq. (31) we get: 

4^^" - 2 (3tanh2e - l) cP^P = Js{0 + 6 - tanh^^) (32) 



for the static part of solution, and: 



- 2 



3tanh^e-l <p'd^=Jd{t,0 
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+6(0(°)^-tanh^e)rf\ (33) 



for the time-dependent part. In the last equation we use the standard nota- 
tion: 

5 Static correction to the polynomial solu- 
tion 



Our goal in this section is to solve Eq. ( [5^ ) and find static part of the 
function 0(^)(t, ^). We adopt the method of solving, which was discussed in 
the previous section. In the first step we put 0^^^ = on the right-hand 
side of Eq. (^) and solve it in a standard way as a linear, inhomogeneous 
equation. Then solution 0^^^ found in this way is inserted to the r.h.s. of Eq. 
( |32[) and in the next step equation with a new source term should be solved. 
Nevertheless, if expression which occur in a new source term as a result of 
this procedure is small we get pretty good approximation in the first step. 
The term in front of (p^p on the r.h.s. of Eq (|3^) is nonzero only on the border 



of the domain wall and quickly tends to zero when | ^ |— > oo. Eventually we 
expect to have a good approximation in the region outside the wall. This 
result seems to be acceptable - inside the domain wall the accuracy can be 
improved by higher degree of the polynomial in Ansatz (|^. Finally we are 
going to solve the linear equation: 

0«"-2[3tanh2e-l]0«=j,(O. (35) 

The solution can be found by the standard Green's function technique. There 
exist two linearly independent solutions of the homogenous part of Eq. (|35|) : 

fi{x) = — 1:2-' 

cosh X 

1 3 3 X 

f2{x) = - sinh2x + - tanhx + - — -5—. (36) 
088 cosh X 

As the Green's function we take: 

G(e,x) = /i(x)/2(oe(e-x) 

-m)f2{xm^-x)-Q{-x)]. (37) 
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The Green's function G{^, x) obeys the condition: 

G{C, = 0, a;) = dla x G {—oo, oo). 
The general solution of Eq. (^) has then the form: 

0«(O=A/i(O + 5/2(0+ / Gi^,x)jsix)dx, 



R 



(38) 



(39) 



where A and B are arbitrary constants. Inserting the Green's function (P7|) 
into the formula (^) we get: 

0«(O =(a - f2{x)js{x)dx] /i(e)+ Mx)js{x)dx] /2(0. 



(40) 

We have to put A = and B = because we are looking for a solution 
generated by source term, not for the homogeneous equation solution. An- 
other reason for keeping i? = is quick (exponential) growth of the function 
/2(x) for I ^ 1^ 00. In the case 5 7^ the solution grows exponen- 

tially and doesn't meet requirements of the perturbative calculus. Combining 
expressions (pHD, (E7\) and (HO) we get finally: 



(1) 



(e>l) 



/■I /27 1 \ 
Jo ^^^^^ [t^^ ~ 3^^^^ ^ 



and 



+ 



-1 



'27 
4 



X X I + X ] dx 



+ x^j dx 

/2(0- 



(41) 

/i(0 
(42) 



Solution for ^ < has opposite sign. 

Integration in the formula (^Tj) can be done numerically, yielding: 



cosh2(e)' 



(43) 



where C ~ —0.731. As one can expect 0s^^(O tends to zero when ^ — »• 00. 
The static solution with the correction 0^°^ + 0^^^ is presented in Fig. 2, 
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where one can compare it with the pure polynomial solution (p^^^ and the 
strict, static solution tanh ^. As it is easy to notice the correction obtained 
above forces the polynomial solution to the well-known strict, static solution. 
It gives us a strong argument, that spliting the static part of polynomial so- 
lution into two parts (one of them is tanh ^ ) is acceptable. This fact will 
be used in the next section. 

Eq. (^2[) can be also solved in a straightforward way using numerical meth- 
ods. Numerical analysis confirms the results obtained above. 



6 Time - dependent correction to the poly- 
nomial solution 

In this section we find time-dependent correction to the polynomial approx- 
imation. We have to solve Eq. (p3D using the iterative method proposed in 
section ^. The second term on the r.h.s. of the Eq. (^) is of the order of 
the small parameter 6 and can be neglected in the first approximation. The 
coefficient in front of 0^^^ in the third term on the r.h.s. is the same as in 



Eq. ([34) discussed in the previous section. Now it is necessary to specify 
more precisely the scenario of evolution of the discussed domain wall. For 
t < the domain wall is static and, as was discussed in the section ^ has the 
form of the strict, static solution tanh ^ (see (|])). In t = the domain wall 
is squeezed by 26 as a result of action of an external force. The evolution at 
later times is governed by the field equation (pSl). 

This scenario helps us to understand the role of the last term on the r.h.s. 
of Eq. (^). In the simplest case we can state that it is small for ^ — > ±00, 
and by neglecting it we obtain good approximation in the first step in this 
region. Nevertheless from the discussion presented above we can conclude 
that we get quite good solution in the whole range of ^. To prove this we 
rewrite the third term on the l.h.s. Thus we have: 



+6 50?) (50?) + 2 50f + 2 tanh ^) (0^^) + 0« j . (44) 

Due to the results of section ^ the static part of the polynomial solution is 
corrected by the function 0^^-*, which has been calculated above. When we 
take some terms of higher order (which can be obtained in the next steps 
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of calculation) we can change (p^^^ ^ on the r.h.s. of Eq. (|33| ) by "improved" 
expression (cpf'^ + (p^^^y . From the discussion presented above and illustrated 
in Fig. 2 we have: 

0(0) + 0(1) -tanh^. (45) 

In the first step we are then going to solve the linear, inhomogeneous equa- 
tion: 



- 2 [3 tanh^ ^ - 1 
This can be rewritten as a wave equation: 

92 



(46) 



(47) 



where operator has the form: 



^-2(3tanh2e-i: 



(48) 



Eq. (47) can be solved by the standard Green's function technique. We cal- 
culate Green's function for opearator using expansion in eigenfunctions. 
We are looking for the Green's function G{C,,t]C,',t'), which fulfill the equa- 
tion: 

G{^,t;^\t') = S{t-m^-a, (49) 



and obeys the condition: 

G{^,t;^',t') = for t<t'. 



(50) 



If the set of eigenfuctions {i'n} is given, we can construct Green's function 
in a standard manner. Detailed description of this procedure is given in 
The Green's function can be written as: 



G{^,t;e,t')=K{^,t;e,t')e{t-t') 
where propagator K{^, t; t) has the form: 



(51) 



An 
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The solution of the inhomogenous equation is given by the formula: 

<l>i\t,0 = - I di' f dt'K{i,t-i'X)j,{t\0- (53) 

J R J —oo 

Because < 0,^) = we can simplyfy (^) as folows: 

4\t,0 = - I di' f dt'K{i,t-i\t')j,{t',0- (54) 

JR Jo 

The problem of calculating time-dependent correction to the polynomial so- 
lution is then reduced to finding the system of eigenfunctions of operator D^. 
This system is well-known (see e.g. )and consists of two discrete eigenvalues 
and continous spectrum with corresponding eigenfunctions: 

for Ao = 0, (55) 
for Ai = 3, (56) 



V3 



2 cosh^ X 
3 sinhx 



2 cosh X 



where: 



Akx 



27r(A;2 + l)(A;2 + 4) 



2 - 



cosh X 



— Siktanhx — 



(57) 



Afc = A;2 + 4, ke R+. 



The eigenfunctions from the continous part of the spectrum correspond to the 
real eigenvalues and can be split into two sets of eigenfunctions, orthogonal 
to each other: 



and: 

^l^\x) 



1 



^27r(A;2 + l)(A;2 + 4) 



2-e 



cosh X 



27r(P + l)(A;2 + 4) 



2-k' 



cosh X 



cos kx + 3k sin kx tanh x 
(58) 

sin kx — 3k cos kx tanh x 
(59) 
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Inserting formulae (^Sj) - (|59|) into (0) - (^ID we get the expression for the 
retarded Green's function: 

Gr{^, t; t') = e{t - 1') {{t - t')MOro{0 

+l=sm[V3it-t')]Mori{a 

+ Er^^ ^"^^^"'^V i-^(Oe(oV (60) 

We get the solution (f)^J\t,^) integrating Gr{^,t;^',t') with the source term 
idii'ii') in the whole range of ^' . As jd(t',^') is odd in variable part of 
Gr even in this variable does not give any contribution to the final solution. 
The Green's function can be rewritten as: 

G,K, = e(i - i') ('-i=sm|v^(i - Ol^Kj^'K') 
+ f (61) 

Finally, inserting Eq. ( pT| ) into (Q) we can write solution (j)l^\t,^). For 
convenience we split it into two parts due to the parts of Green's function 
(|6lD . It reads: 

4\t,o = <p^IUt,o+<i>^IUt,o, (62) 

where and 4>d^l_) are given by formulae: 



(3 cos \ 
(2 — A;^) sin — 3A; cos fcf tanh f — 
^ ^ ^ ^ ^ cosh'e/ 

f (2 - A;2) sin k^' - 3k cos A;^' tanh ^ - ^^S^] ^ (53) 
\ cosh ^' J 

^'(-.'«' - -T^, l^eJ^^mt-oM^sit - (64, 
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Let us give an interpretation of these two parts of the final solution, 
connected with the ipi mode in the system of eigenfunctions of the operator 



.(1) 



D^. The interaction with ipi generates in the time-dependent correction (j)^^^ 
the component, which exponentialy vanishes for large ^ and oscillates. It 
can be treated as a form of excitation of the domain wall given by ipi. The 
crucial role which this component plays in dynamics of our system is more 
transparently visible in the process of collision of kinks (see and |[T7|). 



The component has the form of a wavepacket - the oscillations of the 

domain wall generate radiation. Thus we have the following scenario of the 
discused phenomena: in t = the external force sqeezed the domain wall 
by 26. The wall oscillates around the static configuration. The oscillations 



generate radiation as in the formula (p3|). 
waves going out from the wall: 



It is natural to restrict 



to 




dt' 



1 



(fc2 + l)(P + 4)3/2 

3 



3A;tanh^ [2-k 



cosh ^' 



cosh 



-t')] 

gPtanh^tanh^' 



2 2-k' 



cosh 



(65) 



For 1^ < the result is analogous. 
The solution (p^J^^) has the form of a wavepacket: 



.(1) 
Pd{- 



dkF{k,t,^)e 



(66) 



The waves have the frequency uj{k) = v«?+4 so they satisfy the dispersion 
relation: 

uj\k) -e = A. (67) 
It agrees well with our expectations, bacause for large | ^ | Eq. (|3) reads: 



0. 



(68) 
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To plot the function (jr^ given by the formula (|65|) we calculate it numeri- 
cally for some fixed values of t and ^ G (0, 20), with the step = 0.1. The 
final result is presented in Fig. 3-7. One can easily see the causal character 
of obtained solution. In the region outside the light-cone the field ^^^(L) is 
nearly equal to 0. The fluctuations are caused mainly by computation errors. 
This fact was checked by decreasing the step of integration - the fluctuations 
decrease then too. The field 'P^d}^) small - on the figures it is multiplied by 
factor 10^. This once again confirms the approximation used in our approach. 



7 Backreaction of the radiation 

In the previous section we have shown that the considered, oscillating domain 
wall radiates. Thus, we have the energy flowing out from the domain wall. 
We can then expect dumping of oscillations of the domain wall. In this 
section we construct a simple model of this process. We start from equation 
of energy balans: 

^--m. (69) 

where is the energy per unit area of the domain wall for the polynomial 
solution and F{t) is the energy emitted from unit square in the moment t. 
From Eq. one can calculate Neglecting terms of the second and 

higher order in 5 we get: 

^(0) = Ef)+Ef 5cosfit, (70) 

where Ef'^ ^ 1.58, E^'^ ^ 0.81. Ef^ is the energy of the static, polynomial 
domain wall solution. One can easily compare it with the energy of the strict, 
static solution tanh ^: ii^[tanh ^] ^ 1.33. as we can expect the energy Ef'^ 
is a bit greater than i?[tanh 

In Eq. (^) is time-dependent. It is caused by limitted accuracy of 
the polynomial solution - in the case of a strict solution energy should be 
independent on time. In t = we have Ef"* + 5E^^^ as the energy of our 
configuration. For t > it oscillates as a result of approximation in the 
polynomial solution, nevertheless for the strict solution it should be constant. 
Accordingly we take: 

E(°) = E(°)(0)=Ef . (71) 
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The parameter 6 describes the amphtude of the domain wall oscillations. In 
reahty it is time-dependent: 6 6(t) (and we expect that it decreases as 
a result of radiation emitted from the wall). The energy flux can be easily 
calculated from the energy- momentum tensor for the action (|l|): 

T^^ = 5^05> + r/^^/:. (72) 

The energy flux density is given by the formula: 

ro = 5>9V (73) 



F(t) = -2(9V« ) (74) 



From the continuity equation we can get in a standard manner expression 
for F{ty. 

As ~ (5 we can re-define the field as follows: 

<Pil)=0lU (75) 

Thus, we have: 

F(t) = 52(9°0(^(L)5«0it)) l«=i+^- (76) 
Collecting the formulae (|69D, ([71]) and ([76| ) one gets the following differential 
equation for 6{t): 

m = -6'^ 9«0i^(L)) k=i+. • (77) 

The solution of Eq. ( jn\ ) with the initial condition 5(0) = 6 reads: 

d 

Eventually to find the solution 6{t) one needs the value of the field 4>d^l^) 
on the border of the domain wall, which was calculated numerically in the 
previous section. 

We have performed the numerical computations for three different values of 
the initial oscillation amplitude, namely: 6{0) = 0.01, 0.1, 0.3. The final 
result is presented in Figs. 8, 9 and 10. The dumping of oscillations, as we 
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can expect, is the quickest when we have the greatest amphtude. 
The plots of function 6{t) have characteristic bends at multiphcities of half- 
period of oscillations T = 27^/^7 = 7^ — 0(5^). We can interprete them as 
follows: the energy is most strongly radiated from the wall when the velocity 
of oscillations is the greatest - near the returning points the velocity is small 
and the dumping is slower. This coincides with the half-period of oscillations. 
For small, initial amplitudes of the oscillations the dumping is very weak. In 
the case of 5(0) = 0.01 during the time of approximately 4 oscillation periods 
the amplitude decreases only by 4%. For 5(0) = 0.1 and 0.3 the corresponding 
values are 25% and 75%. These two cases can however be outside the region 
of application of our approximation (it demands 5(0) ^ 1 ). 

8 Remarks 

Let us briefly summarize the main points of our work. We have derived 
the planar domain wall solution with the time- dependent thickness in the 
polynomial approximation. The planar domain wall oscillates with small 
amplitude. We have then proposed an approach enabling us to calculate the 
correction to the pure polynomial solution. We have split this correction 
into two parts. The time-dependent part contains the component which is 
interpreted as the radiation from the oscillating domain wall. 

There are several possible extensions and generalizations of our approach. 
It can be applied without much trouble to the other fleld-theoretical models. 
In our paper we discuss only the case of small oscillations around the static 
solution. As one can easily see in Fig. 1 there exist also solutions oscillating 
with amplitude which is not small. This case should be discussed separately. 
In accordance with our discussion in section |^ we can expect very quick 
dumping of the oscillations. On the other hand, we should remember that 
our method is based on the linear approximation and is reliable for small am- 
plitudes of oscillations only. Thus, the process of creating the kink-antikink 
pair is also possible. 
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Figure captions 

Fig. 1. The configuration space of tlie equation 

Fig. 2. The statical polynomial solution, the improved polynomial solu- 
tion and the strict solution tanh^. 



Fig 


. 3. 


The radiation emitted by the oscillating domain wall for t = 3. 


Fig 


. 4. 


The radiation emitted by the oscillating domain wall for t = Q. 


Fig 


. 5. 


The radiation emitted by the oscillating domain wall for t = 9. 


Fig 


. 6. 


The radiation emitted by the oscillating domain wall for t = 12 


Fig 


. 7. 


The radiation emitted by the oscillating domain wall for t = 15 


Fig 


. 8. 


The dumping of oscillations for the initial amplitude 6 = 0.01. 


Fig 


. 9. 


The dumping of oscillations for the initial amplitude S = 0.1. 


Fig 


. 10 


. The dumping of oscillations for the initial amplitude 6 = 0.3. 
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